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Abstract. We investigate the problem of recovering the initial displacement / 
for a solution n of a linear, isotropic, non-homogeneous elastic wave equation, 
given measurements of u on [0, T] X dQ, v^here f2 C is some bounded domain 
containing the support of /. For the acoustic wave equation, this problem is 
known as thermoacoustic tomography (TAT), and has been well-studied; for 
the elastic wave equation, the situation is somewhat more subtle, and we 
give sufficient conditions on the Lame parameters to ensure that recovery is 
possible. 

1. Introduction 

In this paper, we consider the linear, isotropic elastic wave equation and Cauchy 
initial value problem given by 

'{df + P)u =0 in(0,T)xR3^ 
(1) I u\t=Q = /, 

where u = {ui,U2,U3) is the displacement vector, 

-Pu = V • {fi{x){{Vu) + (Vuf)) + V(A(a;)V • u), 
X and ^, arc the Lame parameters, 

^"^-^ = fe- 

is the Jacobian of u and (Vu)-^ is its transpose. Equivalently, we can write 

3 

-Pu = ^ (All + V(V • u)) + AV(V • It) -f (V • u)VA -t- ^ • (Vu-j + d-ju) e^. 

We will assume that X{x) and fJ.{x) are independent of t; we will also assume that 
they are each positive, in that there is some positive constant ao so that A, > 
ao, and later, we will make further assumptions about bounds and smoothness. 
Throughout, we will assume that our initial-data, /, is compactly supported in a set 
fl, which could, in general, be any bounded domain in R'^, though for simplicity, we 
will take ft ~ -8(0, R) for some R > 0. The objective of thermoacoustic tomography 
in elastic media is to recover /, given the data 

A/ := u|[o,T]xaf2- 

In particular, we will show that, if the Lame parameters satisfy a condition on 
their relative size (briefly, that sup y/JI < 3 inf ^/X + 2/i) and one regarding their 
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gradients (to be discussed later), then / is recoverable via a Neumann series. Our 
method is strongly inspired by the techniques of thermoacoustic tomography for 
scalar wave equations, and we adapt the time-reversal approach used by Stefanov 
and Uhlmann in [19j . We devote the next section to discussing the existing results 
and methods in thermoacoustic tomography, and how they motivate the work in 
this paper. 

2. Background and Thermoacoustic Tomography 

Thermoacoustic tomography (hereafter TAT) is a method of medical imaging 
where the object of interest is exposed to a short electromagnetic (EM) pulse, 
absorbing some of the EM energy in the process. Because cancerous cells will 
absorb more of this energy than healthy cells will, it would be diagnostically useful 
to know the absorption function a{x). To accomplish this, TAT makes use of the 
elastic expansion in the nearby tissue caused by the energy absorption, as well as the 
fact that this initial pressure distribution is roughly proportional to the absorption 
distribution. This initial pressure, in turn, leads to a pressure wave p{t, x) that 
propagates through the object, and is then measured by transducers located on 
an observation surface T surrounding the object for some length of time, the goal 
being to use this data to reconstruct the initial pressure. In much of the literature, 
the intended application is the imaging of cancer in a human breast, though some 
authors have written specifically about imaging the brain. This is somewhat more 
complicated, as the skull introduces a jump discontinuity in the sound speed, but 
under certain assumptions, recovery is still quite possible (see |20j . jl8j). 

Mathematically, we consider the scalar wave equation and Cauchy initial value 
problem 

'{df + A)p =0 in(0,r)xM", 
< p\t=o = /, 

^dtp\t=o =0, 

where A{x,D) ~ — c^(a;)A, (or more generally, A{x,D) = — c^(a;)Ag, for some 
other metric g on M"). The initial pressure function / is typically assumed to be 
compactly supported inside of some bounded domain fl C M" (corresponding to 
the object to be imaged), though in some work, / is taken to be supported in some 
compact set, or even merely f G ioi p > 2n/{n — 1) (see p^). The observation 
surface T is often taken to be dfl (sometimes referred to as complete data), though 
there are also satisfactory results for the case where F is some other set, such as 
a portion of (likewise, incomplete data; see, for instance, [11], [37], [25]). The 
data one then collects is 

Af :={p{t,y):0<t<T,yeT} 

and from this, the goal of TAT is to recover /. 

There are three main recovery methods used in TAT, and their applicability 
depends largely on the assumptions made about the geometry and physical at- 
tributes of the medium (in terms of and c{x)) and the observation surface T. 
We will briefly discuss some of these here; for a detailed account and comparison 
of these techniques, including their relative advantages and disadvantages, see the 
excellent survey papers by Hristova, Kuchment and Nguyen , or Kuchment and 
Kunyansky [12], [13]. 
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In the method of fihered backprojection, one assumes that c is constant (i.e. 
the medium is acoustically homogeneous) and that the observation surface F is a 
sphere of radius R. We can then recover / through integral formulas such as 

/(^) = -^ZTB^ / f7 — Z\ ^^(2^)' 



Stt^R 7r \y - x\ 

where 



are the spherical integrals of /, with dA and duj are the surface measures on the 
respective spheres. For more on the integral geometry approach to TAT, see the 
work of Agranovsky, Berenstein and Kuchment [T], Finch, Patch and Rakesh [S], 
Xu and Wang [26], Finch, Haltmeier and Rakesh [7] and Kunyansky |15) . 

The method of eigenfunction expansion applies to a slightly more general setting, 
in theory allowing c to be variable, and for F to be any closed surface (i.e. F = Oil 
for some fl). We then seek to write / as a Fourier series 

where ipk are eigenfunctions of the operator — c^(a;)A in with Dirichlet boundary 
conditions on d^l, and then find the coefficients fk using integral formulas. In |14| . 
Kunyansky showed, in the case F is a cube and c is constant, that / can be recovered 
fast and precisely. Of course, for a complicated set ft, or a variable speed c{x), the 
eigenfunctions and eigenvalues may not be known, and it is not clear whether this 
method can be effectively implemented. 

The third method (and the method we will eventually use) is known as time- 
reversal, and was first proposed by Finch, Patch and Rakesh in [S], and first im- 
plemented by Burgholzer, Matt, Haltmeier and Paltauf in [3]. Here, c is allowed to 
be variable (i.e. the medium can be non- homogeneous), and the restrictions on the 
geometry of F are far less than those of the other two methods. The key assumption 
is that there is good local energy decay, meaning that for / compactly supported in 
fi, the energy of the solution in $7 decays sufficiently fast as t increases. To better 
illustrate the nature of this requirement, we first discuss an ideal case: the constant 
speed wave equation {df — A)p = in M+ xM.^^. In this setting, Huygen's principle 
would apply, and we would know that for some T, u{t, x) = inside fl for t > T . In 
this case, we could recover / by considering solutions of the initial value problem 

{{dl - A)q =0 in M+ x M", 

q\dn = 9{t,y) 

where g = A/, and then solving the problem in the reverse time direction, since 
(7(0,3;) = f{x) by uniqueness. 

Of course, Huygen's principle does not apply in the general setting, but there 
is a useful analog which still leads to strong results. To ensure the kind of energy 
decay we need, we will assume that the speed is non-trapping, meaning that all 
rays starting in f2 leave in finite time, and that the supremum of these times is 
also finite. Explicitly, for the Hamiltonian H{x,^) — hc^ix) |^|^, we consider the 
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Figure 1 . An example of a trapped geodesic, 
solutions in M^'^ of the system 

with G and 7^ (these solutions are called bicharacteristics, and their 
projections to M" are often called rays). We say that {il, c) is non-trapping if each 
ray leaves fl in finite time, and the supremum of these times is finite as well. For 
a non-trapping (f2,c), we will call the supremum T{fl). 

One may gain additional understanding of the nature of this definition by con- 
sidering an example of a trapping metric: consider c(x) = |x|, with an annulus 
{a; : ri < |x| < r2}. Then, for xq G fl, if is perpendicular to xq, it is straight- 
forward to check that the resulting ray is {x : \x\ = \xq\} (see Figure [T]), which will 
remain inside 51 for all time. 

Because singularities propagate along bicharacteristics, the non-trapping hypoth- 
esis also ensures that any singularities of / will have left fl by T{Q), or put another 
way, that p{t,x) is a smooth function on fl for t > T(fl). This has been discussed 
extensively in the literature; see for instance [5], [9] or [25j . 

Because we can no longer assume that p eventually vanishes inside fl, it would 
be incorrect to think that, for any finite time T, solutions q of ^ will satisfy 
q{0,x) = f{x), at least not exactly. To alleviate this, some authors have replaced 
the boundary condition by q\dn = x(^)ff(^i y)y where x is a smooth cutoff function 
vanishing near T > T{^1), and x = 1 near {—oo,T{n)) (see [10]). In ^9\, Stefanov 
and Uhlmann used a different time-reversal method (this will be discussed thor- 
oughly in the final section of this paper), and showed that, for the acoustic wave 
equation with variable sound speed, it is possible to invert A using a Neumann 
series. Subsequently, in |18| . Qian, Stefanov, Uhlmann and Zhou went on to give a 
specific numerical algorithm for recovery based on this theoretical understanding. 
In this paper, we will follow the method of [19] when possible, assuming that the 
absorption of EM energy and subsequent elastic tissue expansion lead to elastic 
wave propagation instead of acoustic. 
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The differences between scalar equations and systems will introduce some diffi- 
culties, the main difficulty being the unique continuation problem; i.e., determining 
when specifying Cauchy data on a hypersurface 5' C M x is sufficient to uniquely 
determine the solution in a neighborhood of S. For the classical case (where the 
coefficients are smooth), this is Holmgren's theorem (see, for instance, [H])- In the 
case of non-smooth coefficients for a scalar wave equation, there are good results 
due to Tataru (see [21],[22])- For the static Lame system (i.e. the elliptic system 
Pu ~ 0), Lin, Nakamura, Uhlmann and Wang have recently shown (in [TB]) that 
there is a strong unique continuation principle; specifically, if A S L°° and /i G C^'^, 
in M" with n > 2, if a solution is zero in a neighborhood of any point, it is identi- 
cally zero. For the elastic wave equation, however, the situation is somewhat more 
subtle; while some results exist, we will have to make assumptions about the Lame 
parameters to acquire satisfactory results, the topic of the next section. 

3. Unique Continuation and Conditions on the Lame Parameters 

One of the key attributes distinguishing the elastic wave equation from scalar 
wave equations is the presence of two speeds of propagation, the so-called P- and 
S- waves or modes (also known as compression and shear waves or modes). In 
particular, these two speeds ci(x) — ^/2^ + A and C2(x) — correspond to the 
eigenvalues of the principal symbol of P: 

Pix, = (A + 2^i)ix)^e + K^) (lei' I - CC"^) • 

It will be useful to define scalar wave operators corresponding with these speeds, 
and so with oi = l/(2/i + A) and 02 = we define Oaj,j = 1, 2 by 

Da^ a,d^ - A. 

The property of finite speed of propagation for the elasticity system is essentially the 
same as the scalar demonstrated by the following definition and theorem. 

Definition 3.1. We will say that u has finite speed of propagation in (0,r) x 
B{xq, e), with maximum speed c > 0, if for any to G [0, T), w(ioi ■) — ut{tQ, •) = 
in B(xo, e) implies u = a.e. in the cone Uo<s<e/cCs, where Cs ~ {t ^ to + s} x 
B{xo, e — cs). 

Theorem 3.2. Assume that p.,X E C^(M'^), and suppose u £ solves (tlj. Then, 
for any open ball B(xo, e) <Z MP , u has finite propagation speed in (0, T) x B{xo, e), 
with maximum speed c — sup^^^^^^ ■\/2/i -|- A. 

For a proof, we direct the reader to pjTj. With this result in mind, we will 
hereafter assume that there exist constants c~ , c"*" so that 

(3) — sup y2/7+A < 00, 

(4) c' = inf v7^ > 0. 

We will also assume that (17, P) is non-trapping, meaning (recalling the definition 
and discussion for the scalar case given on page|4]) that every ray (i.e. the projec- 
tion of every bicharacteristic starting in VL) leaves £7 in finite time, and that the 
supremum of these times, which we will again call T{Vl), is also finite. 
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The presence of two speeds of propagation causes some difficulty regarding ques- 
tions of unique continuation. For a homogeneous medium (i.e., the Lame parame- 
ters are constant), the system is diagonahzable, and the P- and S- modes are pre- 
served throughout the wave's evohition, effectively reducing this problem to the 
scalar case. In the more general setting however, P-waves may transmit or reflect 
as S-waves (and vice-versa) at an interface, making the question of unique contin- 
uation more subtle. For instance, if the P-wave vanishes on the boundary of a set 
for all time, one cannot necessarily conclude it vanishes on the interior as well; it 
may simply be transmitting as an S-wave instead. The extent to which the two 
modes can be decoupled is useful for understanding the reflection and transmission 
of singularities, and has been studied using the pseudodifferential calculus (see [23] 
and [ST for smooth /i and A, and most recently, [2" for /i, A € C^'^). 

Returning to unique continuation, many satisfactory results have been proven 
via Carleman estimates, and a thorough discussion can be found in the work of 
EUer, Isakov, Nakamura and Tataru, as well as that of Cheng, Isakov, Yamamoto 
and Zhou (see and respectively) and other authors. In the following two 
theorems, u will be a solution of ([T]), f2' is an open domain in R^^, and T' will be 
some positive real number (in practice, we will have fl' containing fi, and T' larger 
than T{il)). Both of these results are from [B]; the former is a slightly simplified 
statement of Corollary 3.5, and the latter is Theorem 5.5. 

Theorem 3.3. (EUer, Isakov, Nakamura, Tataru) 

Let fli = i, a2 = '^"■'^ assume, for some 9 > 0, they both satisfy 

(5) 9'^aj{aj + aj^^'^ KVaj|) < +l/2x- Va^ 
and 

(6) e'^aj < 1 

on [-T',T'] X W, that a.j G C^{{-T' ,T'\ x IF), and that fl' C B{0,6T'). Then, if 
u = and d^u = on {~T',T') x dn' , then u{t,x) = when \x\'^ > O'^t^ . 

Theorem 3.4. Assume that the coefficients n,\ & are time independent. Let 
S he a noncharacteristic surface with respect to both Oa-^ and . Then we have 
unique continuation across S for H^^^ solutions u to 

With the assistance of these theorems, we conclude this section by proving a 
unique continuation result needed later for the reconstruction process. Essentially, 
we seek to answer the following question: 

Suppose a solution u of (P), with u{0,x) — f{x) — outside fi, also vanishes 
outside 17 at some later time T. From this, can we determine if f{x) =0 inside VI 
as well? Because of finite propagation speed, this is certainly not the case for all 
times T > (let T be small, and let the support of / be, for instance, some small 
ball contained in 17), and so we seek to describe sufficient conditions for making 
such a determination. 

Theorem 3.5 (Sufficient conditions for the Lame Parameters). Suppose that ji, A e 
are time independent, and that the maximum and minimum speeds of propagation 
satisfy the inequality 

(7) c+<3c-. 
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that there exist 9, T, and e > so that 

(8) 

(9) T> ^(^ + ^) 

and that (ai,a2,T, 0) satisfy the gradient condition ^ on 

3T 3T 



T 

X E ( 0,i?+ -c^ 



Assume also that the surface S — {{t,x) : \x\^ =eH^^ is non- characteristic for 

□ai , ■ Then, for solutions u of ^ (with f compactly supported inQ ^ B(0, R) ), 
if u{T, x) = Q for x ^ fl, we have f{x) = 0. 

Proof. By assumption, we know 

u{T, x) = 0, for x ^ n, 
and since / is compactly supported in fl, we know 

u(0,x)=0, ioi-x<^n 
as well. Thus, by finite speed of propagation, we then have both 
u{t, x) = when |a;| - i? > c+ \T - t\ 

and 

u{t, x) ^0 when |a;| - i? > c+ \t\ . 
Combining these observations shows (see figure [5]) 

T 

u{t, x)=Q when \x\-R> -c+, -~T/2 < t < 3T/2. 

Next, by time-reversal, u extends to an even function of t, so, in fact, we have 

T 

u{t, x) = when - i? > -c+, -3T/2 < t < 3T/2. 

Now, because 9 < c^ , we have 9^aj < 1 for j = 1, 2, and by hypothesis, the gradient 

condition is satisfied on B (0,R+ -jc+ + e) x [— ^] ■ Because 9 > c+/3 and 

T >2{R + e)/{39 - c+), we check that 

T T T 3T 

(10) + _e+ + , < -(30 - c+) + -c+ < —9 

so that we can apply Theorem [331 with 17' = 5(0, i? + f c+ + e) and T' = 
Thus, u{t,x) = whenever \xf > 9^t^, showing that u{0,x) = 0, except possibly 
at the origin. Theorem 13.41 allows us to extend this solution uniquely, and thus 
u{0, 0) = as well, showing / = 0. □ 

Before we continue, we make a few remarks about the conditions placed on the 
Lame parameters in this theorem. First, the condition c+ < 3c^ is needed to 
guarantee we can choose 9 so that ^c^ < 9 < c~; we need 6* < to ensure 
condition ([5]) of Theorem 13.31 is met, and ^c^ < 9 is needed to ensure that the 
first inequality in ((TU]) is valid. This condition essentially reads, "the speed of 
the S-wave must be more than a third that of the P-wave." For many materials, 
this is a realistic assumption (for instance, this is true of most earth materials). 
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Figure 2. The solid gray region is B{0,R) x T); the dashed 
region where w = by finite speed of propagation and time- 
reversal. 



though as McLaughlin and Yoon note in p2], in biological tissue, this assumption 
may not be reasonable. In spite of this, because the observation surface is not a 
true discontinuity in the medium, it may be reasonable to expect that there is no 
exchange between the P- and S-modes, allowing us to work with a stronger form 
of unique continuation and drop the assumptions on the wave speeds. We will not 
pursue this issue any further in the present work, however. 

The gradient condition is somewhat more subtle, though making some additional 
hypotheses could help to simplify the situation. For instance, if we were able to 
assume that the medium is homogeneous outside some neighborhood of (i.e., A 
and /i are constant, so that Voj = 0), then the gradient condition simply reads 
6 < c~ outside that neighborhood. Alternately, if we were to make assumptions 
on the size of |Vaj|, say |Vaj| < Saj for some ^ > 0, then the condition could be 
reduced to 

9^{a,+afs\t\)+^-\x\<l. 

Again, at this point in time, we will leave such additional assumptions for future 
work, and continue with the hypotheses in the form stated in the theorem. 

4. Energy of Initial Data and Solutions 

As we have previously discussed, the success of the reconstruction process de- 
pends on good local energy decay, ensured by the non-trapping condition. Before 
we move on to the reconstruction process, we will find it useful to specifically define 
two kinds of energy spaces associated with ([T]); first for the initial data, and then 
for solutions. In what follows, C/ is a domain in R'^ (in practice, we will have U — 
or U = R'^). We will begin by defining an inner product for the first space. 

Definition 4.1. For functions f,g, let 

if,9)HoiU)= I A(V-/)(V-.g)-|-Mtr((V/)(Vg) + (V/)^(V.g)) dx, 
Ju 
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where tr indicates the trace. Specifically, we have 

It is easy to see that (•, ■)ho(U) symmetric and bilinear. Furthermore, 
Lemma 4.2. For functions f,g G C^{U), we have 

{f,g)Hn{U) = {Pf,g)L'-'{U) = {f,Pg)L^U)- 

Proof. Using the identities 

(V • (m(V/))) • g = V • ((m(V/))5) - tr(V/)(V5), 

(V • (a*(V/)^)) • .9 = V • ((M(V/)^)g) - ^i tr(V/)^(Vg), 

and 

(V(A(V • /))) • .9 = V • (A(V • - A(V • /)(V • 5), 

we have 

{Pf,9)LHu)^ I A(V-/)(V-.9)+Mtr((V/)(V.g) + (V/)^(Vg)) dx 
Ju 

- / V • mi^f) + iVff)9)) + V • (A(V • /).9) dx 
Ju 

Since g £ C^{U), the second integral is zero by the divergence theorem, and we 
have {Pf,g)L^u) = if,g)Ho{u)- Similarly, using / £ C^iU) will show {f,g)Ho{U) = 
{f,P9)LHu)- □ 

Next, we define H£,{U) to be the completion of C^{U) under the norm 

(this quantity is essentially the total elastic energy of / in U, with a factor of 1/2 
omitted for convenience), and then define the energy space 

n{u) = Hd{u) ® L\u) 

with the norm 

\\{fiJ2)Kiu) = \\fi\\HMu) + \\Mmu)- 

Finally, we define the total energy for a function u{t, x) at time t as 

Euiu,t) = Oll^^fcj) + \\utitr)\\l^u) ' 

i.e., the total energy is the sum of the elastic and kinetic energies of u at time t. 
Before moving on, we note that, as a consequence of these definitions, if u solves 
(HI with / e H]:){n), then we have 

Em^{u,0) = ||/||^^(R3) = ll.fllHo(f2) =En{u,0). 
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5. Reconstruction 

Now that we have ah the necessary prehminaries in place, we are in the position 
to investigate the reconstruction process. As suggested by our earher discussion 
of techniques in TAT, one possible approach to reconstruction would be the use of 
time-reversal for solutions vq of 

'{df + P)vo =0 in(0,T)xf], 

Vo\t=T = 0, 

,dtVo\t=T =0, 

where h = A/, and then attempt to define an inverse by 

Aoh -.^ vo{0, ■) inf^, 

with the hope that AqH — AgAf approximates /. The problem with this method 
is that h may not vanish on {t = T} x 90, causing the boundary conditions to be 
incompatible. 

To correct this, we will need to introduce an error term, and modify our approach 
accordingly. Given h (again, eventually h — A/), define (j) and v by 

Pcj) = 0, (j)\an = h{T, •), 



(12) 



'{df + P)v =0 in(0,T)xf7, 

^^l[o,T]xan = 

v\t=T = <t), 

dtv\t=T = 0, 



and note that the boundary data is now compatible to first order, in that (p = h on 
{t — T} X dil. Now, we define the pseudo- inverse A by 

Ah v{Q, •) in n. 

Our goal now is to show that A maps the range of A to i?£)(ri), and that the error 
Hoin) small compared to \\f\\Ho{n)' particular, we will show that 
the operator K — {I — AK) is a contraction on Hd{^). The following theorem is 
our main result. Recall that we say (fi, P) is non-trapping if every ray starting in 
Q. exits in finite time, and the supremum of these times (which we denote by T(f7)) 
is finite as well (see page |4] for more discussion). 

Theorem 5.1. Suppose that f is compactly supported in fl, that {^,P) is non- 
trapping and satisfies the hypotheses of Theorem \3.5[ with T > T{Q). Then AA = 
I — K , where K is compact in H]j{^) and ||^||/fj-,(f2) < 1- ^'^ particular, I — K is 
invertible on Hjjiyi), and we have the following Neumann expansion for f : 

oo 

/ = ^ K'^Ah, h := A/. 

m=0 
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Proof. Let w solve 

'{df + P)w =0 in(0,r)xl7, 

w\[o,T]xan =0, 

w\t=T = u\t=T - 0, 

dtw\t=T =dtu\t=T- 



(13) 



Let w be the solution of ([12]) with h — A/. Then v + w will solve same initial 
boundary value problem as u (with initial conditions given at f = T), and thus 
u = V + w. For t = 0, we then have 

f^AAf + w{0,-), 

and so, 

Kf = w{0,-). 

For convenience of notation, let — u{T, •) and uj^ — ut{T, •). Because = (j) 
on dV,, and because — 0, 

(u^ - 0, (p)Hn{n) = (u^ - 0, P4>)L^{n) = 0. 

Thus, 



I T _ ,\\'^ _\\ T||2 _ <^ II T||2 



showing that 



Because w = on dft, we have {w, w)ijj3(si) = (^'fi', w) L^{n) so that 

= {Pw',w)L2(n) + {Pw,w')L2(n) + (w", u;')l2(o) + (w',w")L2(n) 

= (Pw', + (w',w")l2(o) 

= (w',Pw)i2(o) + (■u;',w")l2(o) = 

(note that this is just a reflection of the fact that the Dirichlet boundary conditions 
imposed on w are energy-conserving). Therefore, 

En{w,0) = Eniw.T) < En{u,T) < E^3{u,T) = En{u,0) = ||/||?j^(a) , 
showing 

We next seek to show that this inequality is, in fact, strict. So, suppose that 
there is some / so that \\K f\\f^^(^^^~^ = This implies that all the above 

inequalities are actually equalities as well, and in particular, Eq(u, T) — E^3 (it, T), 
so that 

u{T, x) = 0, for X (f. n. 

By Theorem 13. 5| we then have f — 0, and thus < ll./ll_ffD(o) ^^^r aU 

non-zero / G ^^^(fi). 

Now, we will show that K is compact. We have that u(T, •) and ut{T, ■) are 
smooth as functions on because T > T{il) (all singularities beginning in ft will 
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have left), and so, as linear operators acting on /, they have smooth Schwarz 
kernels. This will also imply cj) is smooth, by elliptic regularity. From this, we can 
see that the map 

has a smooth kernel, and is therefore compact. Next, we note that the solution 
operator of ([T2]) from t = T to t = (i.e., the map {u'^ — (j), uj) ^ (w(0, •), Wt(0, •)) 
is bounded in "HiU) (unitary, actually, because E^{w,T) — Efi(w,0)). Thus, we 
conclude that the map K : H£){^1) — H]j(n),f i— w{0,-) is compact, as the 
composition of a compact and a bounded map, finally allowing us to conclude 
ll-^llffo(si) 1' ^^'^ proof is complete. □ 

It is additionally worth noting that, in addition to uniqueness, the proof of this 
theorem also gives insight into the stability of this problem. In particular, we have 

providing a bound on in terms of the local energy decay. Furthermore, 

because / — AAf = Kf, this also describes the error in the reconstruction if we 
only use the first term K'^Ah ~ AAf of the Neumann series. Additionally, if T is 
chosen so that the hypotheses of Theorem 13.51 are satisfied, but T < T($7), then 
we can still conclude that \\Kf\\^^^^.^ < for all non-zero /, but we can 

no longer be sure K is compact, and thus have no reason to expect ||^||/f^(f2) 
strictly less than 1. Finally, if T < then one can find a non-zero / so that 
ll^/ll_f/D(n) = ll/llHD(n) (by choosing an / supported in a small enough ball, finite 
speed of propagation will imply u{T,x) — outside of f2). 

6. Conclusion 

While we have suggested a method of reconstruction for thermoacoustic tomogra- 
phy in elastic media, and demonstrated the uniqueness and stability of this method, 
there is still much work to be done. 

Indeed, we have only explored the case of complete data, and in practical appli- 
cations, such as breast imaging and geophysical imaging, one can only expect to be 
able to measure data on a portion of the boundary. Thus, examining the problem 
of incomplete data will be critical to the practical relevance of the elastic approach, 
and should be explored in future work. As previously mentioned, various authors 
(again, see [H], [57], [5S]) have investigated this problem for the acoustic wave 
equation, and produced satisfactory results; perhaps by adapting these methods, 
good results can be obtained for setting of elastic media as well. 

Furthermore, we have worked with the restriction that < c^, which, as 

mentioned, is a reasonable hypothesis for earth materials, but not necessarily for 
biological tissue. This assumption may be stricter than actually needed, and so it 
would be nice to find an alternative to using [5| in the uniqueness step. At this 
point, it is not clear how to accomplish this, but perhaps by adapting the methods 
(rather than the results) of Eller, Isakov, Nakamura and Tataru some progress can 
be made. Indeed, in their work the coefficients aj are generally allowed to vary in 
time and be of relatively limited smoothness (e.g., C^); possibly by relaxing these 
assumptions, a better unique continuation result can be obtained. 
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Finally, developing computer simulations and attaining numerical data based 
on the method proposed in this paper are an obvious next step for this approach 
to mature from theory to application. As we discussed earlier, this has already 
been done for the acoustic case (in |18) and others), and hopefully computational 
methods can be developed for elastic media in the near future. 
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